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ABSTRACT 


A feature of great interest in the late Eocene-early Oligocene volcanic province of the Sierra Madre 
del Sur is an elliptical NW-SE oriented dome structure (52 x 30 km) in the Coxcatlan-Tilzapotla region. 
The elliptical structure encompasses the Tilzapotla collapse caldera, rhyolitic domes, large volumes of 
ignimbrites, as well as the Buenavista intrusive body, and the Coxcatlan and Chautle plutons located west 
and east of the structural margin of the caldera, respectively. Previous geochronological studies carried out 
on the silicic and intermediate magmatic rocks intercalated with pre-caldera tilted terrestrial sediments, 
as well as the occurrence of an angular unconformity separating fine-grained and coarse-grained beds, 
constrain uplift in the dome area to the late Eocene (~38-34 Ma). This suggests that doming was related 
to emplacement of magmas into the crust prior to collapse of the Tilzapotla caldera at 34.3 Ma. 

New LA-ICPMS zircon U-Pb and “Ar-*°Ar mineral ages for the key magmatic units combined with 
field observations, petrographic studies and geochemical analyses reveal the geochronology of magma 
emplacement and development of the volcano-plutonic complex, and its connection with the dome. 
Zircon *°Pb?*U ages of the Coxcatlan pluton ranging from 39.5 + 0.6 to 35.2 + 0.2 Ma suggest that 
the Coxcatlan pluton was constructed gradually over ~4 m.y. and overlapped in time with the beginning 
of the development of the Tilzapotla ignimbrite (36.6 + 0.4 to 33.3 + 0.5 Ma). Significant overlap of the 
population density distributions is found indicating temporal continuity of the zircon crystallization 
record. Zircons from the mafic Chautle pluton (34.37 + 0.26 Ma) yield the same weighted mean age within 
uncertainty as zircon from the Tilzapotla ignimbrite. These data are best reconciled with a continuous 
but incremental assembly of a volcano-plutonic system that climaxed with the caldera collapse and 
the eruption of the Tilzapotla ignimbrite. The evidence of mafic inputs throughout the history indicates 
that recharge provided the thermal engine for the prolonged history of the system. Mutual intrusion 
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relationships observed between mafic enclaves and granodiorite facies indicate magma mingling within 
the magma chamber when both host and enclaves were not completely crystallized. 

Development of the Coxcatlan-Tilzapotla system resulted in updoming of the surface, extending 
the “roof” above the melt-rich part of the system in the southeast. The spatial and temporal relationship 
between the crystal-rich Tilzapotla ignimbrite and the less differentiated Chautle intrusion suggests 
that mafic magma injection might have disrupted crystal mushes in the chamber, and destabilized the 
“perched” magma reservoir and thereby triggering the eruption in the southeastern part of the dome. 


Key words: volcanic-plutonic connection, ignimbrite, Coxcatlan pluton, Tilzapotla caldera, southern 
Mexico. 


RESUMEN 


Un rasgo de mucho interés en la provincia volcanica del Eoceno tardio-Oligoceno temprano de 
la Sierra Madre del Sur es la estructura domica de forma eliptica, orientada al NW-SE (52 x 30 km) y 
ubicada en el area de Coxcatlan-Tilzapotla. La estructura eliptica contiene la caldera de Tilzapotla, domos 
rioliticos, volumenes importantes de unidades ignimbriticas, asi como el intrusivo de Buenavista y los 
plutones de Coxcatlan y Chautle, ubicado estos ultimos al occidente y al oriente del anillo de la caldera, 
respectivamente. Los estudios geocronologicos realizados previamente en las rocas magmaticas silicicas 
y las intermedias intercaladas en sedimentos continentales anteriores a la caldera, asi como la presencia 
de una discordancia angular entre capas de grano fino y grueso, constrinen el tiempo del levantamiento 
del domo al Eoceno tardio (~38—34 Ma). Lo anterior sugiere que el proceso de formacion del domo 
esta relacionado al emplazamiento de magmas en la corteza previo al colapso de la caldera (34.3 Ma). 

Nuevas edades de U-Pb por LA-ICPMS en zircones y edades de “Ar-**Ar en minerales de unidades 
magmaticas clave, combinadas con observaciones de campo, estudios petrograficos y andlisis geoqui- 
micos ponen de manifiesto la geocronologia del emplazamiento magmatico y el desarrollo del complejo 
volcanico-plutonico, asi como su conexion con el domo. El rango de edades *°Pb/*U en circones del 
pluton de Coxcatlan de 39.5 + 0.6 a 35.2 + 0.2 Ma sugiere que el pluton de Coxcatlan se construyo gra- 
dualmente durante ~4 m.a. y se traslapa con la historia magmatica de la ignimbrita Tilzapotla (36.6 + 
0.4 a 33.3 £0.5 Ma). Se observa un traslape significativo en la distribucion de edades de las poblaciones 
que indica una continuidad temporal en el registro de la cristalizacion de los circones. Los fechamientos 
de circones del pluton mafico de Chautle arroja una edad igual (34.37 + 0.26 Ma), dentro del error, a 
las edades principales de la ignimbrita Tilzapotla. Estos datos se ajustan mejor con un proceso continuo 
pero creciente de construccion de un sistema volcanico-plutonico que alcanzo su climax con el colapso 
de la caldera y la erupcion de la ignimbrita Tilzapotla. La evidencia de la inyeccion de magmas mdficos 
durante este proceso indica que la recarga de magmas proporciono la energia térmica durante la historia 
prolongada del sistema magmatico. Las relaciones mutuas de intrusion observadas entre los enclaves 
maficos y la facies de granodiorita indican mezcla inhomogénea (mingling) dentro de la camara cuando 
ambas fases no habian alcanzado Ia cristalizacion completa. 

El desarrollo del sistema Coxcatlan-Tilzapotla result6 en la formacion del domo, extendiendo 
el techo arriba de la zona rica en fundido hacia el sureste. La relacion espacial y temporal entre la 
ignimbrita Tilzapotla, que es rica en cristales, y la intrusion menos diferenciada de Chautle sugiere que 
la inyeccion de magmas mdficos en la camara podria haber perturbado las acumulaciones de cristales 
(mushes) y desestabilizo el reservorio magmatico en la parte superior de la corteza causando la erupcion 
en la parte sureste del domo. 


Palabras clave: conexion volcanico-plutonica, ignimbrita, pluton de Coxcatlan, caldera de Tilzapotla, 
sur de México. 


INTRODUCTION 


There has been considerable debate in recent years 
about the relationship of large caldera systems and batho- 
liths. The debate has centered on the largest volcanic 
systems, commonly called “supervolcanic”. It is generally 
agreed that large silicic eruptions can be interpreted in terms 
of “snapshots” that record conditions in the magma chamber, 


whereas plutons provide a longer record of magmatic evo- 
lution (Bachmann et al., 2007b; Lipman, 2007). However, 
some view these as separate systems where batholiths 
represent slow accumulations of silicic magma that never 
attained eruptibility, while supervolcanic systems represent 
rapid accumulations of silicic magma that have no plutonic 
record (Glazner et al., 2004; Coleman et al., 2004, Zimmerer 
and McIntosh, 2012; Tappa et al., 2011). The alternative 
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view most recently articulated by the likes of Bachmann et 
al. (2007b), de Silva and Gosnold (2007), Lipman (2007), 
and Miller (2008) is that large silicic eruptions associated 
with collapse calderas are intimately related to contempora- 
neous, consanguineous, co-genetic intrusives, that represent 
“non-erupted” dominant volumes of the magmatic systems 
(e.g., Smith, 1979; Lipman, 1984; Elston, 1984). In this 
view, individual calderas record the development of a pluton 
and different cycles of activity from that system would mark 
the progressive or incremental development of that pluton. 
Collections of spatially and temporally-associated calderas 
would record the development of a composite batholith at 
depth. Herein we provide insight into this debate through a 
study of the Coxcatlan-Tilzapotla region of the Sierra Madre 
del Sur in Mexico. 

The north-central Sierra Madre del Sur (SMS) has 
been recognized as a site of broad exposures of volcanic and 


100°00" 


& QValle de Bravo 
he 
é 


19°00" 


a *. 


Nanchititla — 


La/ ¢ 
a Qa Muiieca ¢ 


Goleta 


18°30' 


San Vicente 


Gulf of 
= Mexico 


1 Mexican Volcanic Belt 
Sierra Madre del Sur volcanic rocks 
[lj Sierra Madre Occidental volcanic rocks 


dome structure 


ib 


Taxco’ 
8 


aS 


plutonic rocks in the central part of the late Eocene-early 
Oligocene silicic volcanic province in southern Mexico 
(Moran-Zenteno et al., 1999). Of particular interest is 
the Coxcatlan-Tilzapotla area, given its location within 
a NW-SE oriented elliptical dome structure (52 x30 km) 
that encloses the Tilzapotla caldera, medium-grained 
plutons, hypabyssal intrusive bodies, rhyolitic domes, 
and voluminous ignimbrites (Figures | and 2). Age and 
compositional similarities in these broad exposures of 
volcanic and plutonic rocks suggest a genetic relationship 
between them, and provide an excellent opportunity for 
studying the evolution of the magmatic system within the 
elliptical dome structure and understanding the genetic 
connection between the volcanic and plutonic realms. 
Based on the geochronologic relationships, compo- 
sitional variations and geochemical characteristics of the 
voluminous late Eocene-early Oligocene volcanic succes- 
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Figure 1. Silicic volcanic centers of the northern Sierra Madre del Sur. Indicated are the elliptical dome structure, Tilzapotla caldera, and near 
N-trending Laramide structures. Inset indicates the location within Mexico and the large volcanic provinces: Mexican Volcanic Belt (Miocene- 
Recent), Sierra Madre del Sur (late Eocene-early Oligocene), and Sierra Madre Occidental (Oligocene-early Miocene). Modified from Moran- 


Zenteno et al. (2004). 
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Figure 2. Geological map of the Coxcatlan-Tilzapotla area. Indicated are the Coxcatlan and Chautle plutons, the Buenavista intrusion, the Tilzapotla 
caldera ring structure, and related tectonic features. 


sions, Mori et al. (2012) interpreted the origin of the silicic de Silva et al. (2006) and de Silva and Gosnold (2007) and 
centers in terms of the gradual thermomechanical maturation are thus interpreted as the sequence of events and conditions 
of the crust due to a sustained mantle power input during the that trigger the formation of ignimbrite flare-ups. 

Eocene. This group of silicic centers displays volume-time Regional geological and geochronological studies 
patterns similar to those documented and summarized by were carried out previously by our group near the ellipti- 
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cal dome structure in the silicic magmatic rocks, the tilted 
pre-volcanic continental sediments and intercalated volcanic 
rocks, as well as the underlying deformed Cretaceous ma- 
rine sedimentary rocks. Results suggested uplifting in this 
region during the late Eocene (~38—34 Ma) and a doming 
process related to the emplacement of magma into the up- 
per crust prior to extrusion of the Tilzapotla ignimbrite at 
34.3 Ma, which is associated with the caldera collapse 
(Moran-Zenteno et al., 2004; Moran-Zenteno et al., 2007). 
Dome construction was initially interpreted by Moran- 
Zenteno et al. (2004) as the expression of the tumescence 
stage in the development of the Tilzapotla caldera; however, 
results presented here indicate a more complex scenario. 

Earlier studies in the Coxcatlan-Tilzapotla sector 
focused mainly on the volcanic rocks and pre-volcanic 
continental sedimentary rocks, whereas plutonic and hyp- 
abyssal intrusive rocks exposed along and outside of the 
Tilzapotla caldera ring, the structural margin of the caldera, 
have received less attention (Figure 2). 

The present study is based on new LA-ICP-MS 
zircon U-Pb and Ar-Ar mineral ages for several of the 
different magmatic facies within the elliptical dome struc- 
ture, with emphasis on the plutonic and hypabyssal rocks 
(Coxcatlan and Chautle plutons, Buenavista intrusion) in 
order to identify intrusive rocks that are contemporane- 
ous with the Tilzapotla ignimbrite and caldera collapse, 
and gain insight on plutonic-ignimbrite connections. We 
include description of the petrography of different mag- 
matic facies of the Coxcatlan pluton, including magmatic 
enclaves, the Tilzapotla ignimbrite and related units. In 
addition we document the temporal evolution of magma 
emplacement within the elliptical dome. Particular atten- 
tion was focused on the development of the elliptical dome 
structure and its possible relationship to the progressive 
construction of a batholith at depth under the Coxcatlan- 
Tilzapotla region. 


GEOLOGICAL SETTING AND STRATIGRAPHY 
OF THE COXCATLAN-TILZAPOTLA AREA 


The Tilzapotla caldera area is part of the subduction 
related magmatic record in the SMS that extends over a 
broad region in southern Mexico and became extinct be- 
tween the early Oligocene and early Miocene (Schaaf et 
al., 1995; Moran-Zenteno et al., 1999). Magmatic activity 
in the Coxcatlan-Tilzapotla sector formed part of a major 
episode of late Eocene to early Oligocene explosive silicic 
activity preceded by pulses of andesitic and mafic volcanism 
in the north-central SMS. Silicic volcanism produced vo- 
luminous ignimbrites, collapse caldera structures, rhyolitic 
domes, pyroclastic dikes, hypabyssal intrusive bodies and 
intermediate to silicic lavas between Nanchititla and Huautla 
(Figure 1). General rejuvenation in the region led to the 
partial exhumation of most of the Paleogene silicic centers 
and the exposure of feeding conduits and plutonic counter- 


parts. After the Paleocene volcanism died out at ~28 Ma, arc 
volcanism reinitiated at 20 Ma in a more northern position 
in the Trans Mexican Volcanic Belt (Ferrari et al., 2012). 

Silicic volcanism of the Coxcatlan-Tilzapotla area 
covers folded Cretaceous marine and sedimentary units 
that include calcareous and clastic beds of the Morelos/ 
Huitzuco and Mexcala formations (Fries, 1960); it also 
overlies Eocene terrestrial deposits with intermediate and 
mafic dikes and lavas intercalated in the sequence (Balsas 
Group) (Fries, 1960; Moran-Zenteno et al., 2007) (Figure 2). 

The structural dome associated with the Tilzapotla- 
Coxcatlan magmatic system is a NW-SE oriented elliptical 
structure with a long axis of 52 km and a short axis of 30 
km (Figure 1). The 30x20 km collapse caldera structure 
is contained in the southeastern part of the dome, whereas 
the Coxcatlan pluton is located in the northwestern part. 
Compared with pre-volcanic ~N-S folds characteristic of 
the Late Cretaceous - Paleocene shortening in the region 
(Laramide), the Coxcatlan-Tilzapotla dome clearly displays 
a different orientation of its long axis. The effects of this 
structural interference are expressed by the northward and 
southward plunging of the small-scale Laramide folds on 
the northern and southern flanks of the dome, respectively 
(Figures | and 3). Timing of the dome formation has been 
constrained by an angular unconformity and a significant 
facies change between two Eocene terrestrial successions 
exposed outside the northwestern edge of the structural 
dome. The lower and upper successions are separated by 
ignimbrites dated by K-Ar in sanidine at 38.5 + 0.7 Ma 
and 35.1 + 0.4 Ma (Moran-Zenteno et al., 2007). Tilting 
of the lower succession, as well as a coarsening-upward 
tendency observed in the upper part of the lower succes- 
sion, is indicative of rejuvenation of the landscape due to 
uplift of the Coxcatlan-Tilzapotla dome structure. The oc- 
currence of ignimbrites as old as 38.5 Ma between the two 
successions and tilting of the upper succession away from 
the dome structure indicate that uplift had already initiated 
at least 38 m.y. ago and continued until ~34 Ma, when the 
Tilzapotla caldera collapsed. 


Volcanic stratigraphy 


Ignimbrites associated with caldera collapse, post- 
caldera ignimbrites and rhyolitic domes in the Tilzapotla 
area were erupted from 35.5 to ~33 Ma (Moran-Zenteno et 
al., 2004). The pre-silicic volcanic record is represented by 
low volume andesitic to mafic lavas and dikes intercalated 
with and intruding the terrestrial deposits of the Balsas 
Group (Fries, 1960; de Cserna and Fries, 1981), which 
include the Tepetlapa and Huanjintlan formations cropping 
out northwest of the dome structure (Moran-Zenteno et al., 
2007). Intercalated with the Balsas Group and separating the 
Tepetlapa and Huanjintlan formations are ignimbrite flows 
dated at 38.5 + 0.7 and 35.1 + 0.4 Ma; the earlier ignimbrite 
has the same age within uncertainty as basal sanidine-rich 
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Figure 3. Digital elevation model of the elliptical dome structure, with cartoon showing the collapse caldera in the southeast sector of the dome. Stereonet 
diagram plots small scale Laramide fold axes north (white circles) and south (black circles) of the dome. Modified from Moran-Zenteno et al. (2004, 2007). 


ignimbrites of the neighboring Taxco volcanic center, sug- 
gesting that the two flows are related (Moran-Zenteno et al., 
2007). The upper dacitic to andesitic lava flows and dikes 
in the caldera area span from 33 to 32 Ma and are related 
to resurgence that produced uplift of the central block of 
the caldera infill (Figure 2) (Moran-Zenteno et al., 2007). 

The Tilzapotla ignimbrite (Figure 4a), which is inter- 
preted as the pyroclastic unit related to the collapse, varies in 
exposed thickness from at least 600 m in the caldera interior 
to 50-150 m of preserved ignimbrite sections outside the 
caldera (Moran-Zenteno et al., 2004). The intra-caldera 
ignimbrite succession is formed by at least five flows that 
display differences in the relative abundances of lithic and 
pumice components. Phenocryst abundance in pumice frag- 
ments varies typically from 25% to 35%, although in basal 
units of the northern extra-caldera facies the abundance is 
less than 5%. The main phenocryst phases in order of abun- 
dance are quartz, plagioclase, sanidine and biotite. Traces 
of altered hornblende are also present. The groundmass of 
the ignimbrites is generally devitrified and is composed of 
microcrystalline quartz and feldspar. The caldera ring on 
the north and south sides coincides with NW-SE trending 
faults that record post-collapse lateral displacement (Figure 
2). The ignimbrite facies near the ring has relatively high 
abundance of limestone lithics and some sections display 
collapse megabreccia. We therefore interpret the “ring” as 
the structural margin of the caldera. 

Post-collapse silicic units include the Rodarte and 
Atopula ignimbrites that overlie debris flows and con- 


glomerate deposits derived from the Tilzapotla ignimbrite. 
The Rodarte ignimbrite contains large crystal-poor pumice 
fragments, whereas the Atopula unit displays crystal-rich 
(30-40%) pumice fragments. Relative phenocryst phase 
abundances are similar to the Tilzapotla ignimbrite, with 
quartz dominating over plagioclase and sanidine. In the cen- 
tral area of the caldera these units are overlain by the Gallego 
formation, a ~50 m succession of vitrophyres interpreted 
as having been formed from densely welded crystal poor 
ash flows. A series of rhyolitic domes with characteristic 
flow structure are distributed along the western segment 
of the caldera ring. 

The post-ignimbrite volcanic record is represented 
by a succession of hornblende dacite to pyroxene-bearing 
andesitic lava flows and feeding dikes, up to 200 m thick, 
that are distributed in the central and northwestern sectors 
of the caldera. The succession includes the dacitic to andes- 
itic El Salto formation and the andesitic Rodeo formation 
(Moran-Zenteno et al., 2004). Differential erosion of the 
caldera and uplift related to the resurgence led to inverted 
relief, so that the top of the caldera infill ignimbrites and 
lavas are now at higher elevations than the surrounding 
Mesozoic calcareous rocks (Moran-Zenteno et al., 2004). 


Plutonic record 


The intrusive record is represented by three plutons 
located outside of, but relatively close to, the caldera struc- 
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Figure 4. Field photographs from the Coxcatlan-Tilzapotla area. a) Tilzapotla ignimbrite containing slightly flattened pumice; b) Mafic magmatic enclave 
(MME) hosted in the Coxcatlan pluton. Rounded form suggests plastic deformation within a fluid host. Double magmatic enclave (darker enclave within 
lighter enclave) suggests mingling of magmas of different compositions. Note reaction rim in contact zone between the MME and Coxcatlan pluton; 
c) Felsic magmatic enclave in the Coxcatlan pluton; d) Lobate form of MME in the light gray colored Coxcatlan granodiorite; e) Aligned belt of MME 
observed in the southern part of the Coxcatlan pluton; f) Double magmatic enclave in the Coxcatlan pluton; g) Felsic and mafic enclaves side by side; 
h) Magmatic enclave swarm in the southern part of the Coxcatlan pluton. 


254 Martiny et al. 


tural margin and within the structural dome containing the 
Tilzapotla caldera. Their position within the dome and the 
spatial and temporal proximity to the Tilzapotla caldera 
strongly suggest that they belong to the same magmatic 
system that produced the caldera. The largest exposure, 
located west of the structural margin, corresponds to the 
Coxcatlan pluton (Figure 2), which represents a shallow 
apophysis of the deeper silicic reservoir in the uncol- 
lapsed part of the dome that has been exposed by drainage 
incision. The other two plutons are represented by small 
exposures located along the structural margin southwest 
of Buenavista (Buenavista intrusion) and east of the struc- 
tural margin (Chautle pluton). The characteristics of the 
plutonic exposures in the vicinity of the Tilzapotla caldera 
were poorly known and only brief mention has been made 
of the intrusive rocks in a local publication and a technical 
report (Ontiveros-Tarango, 1973; Carrizales-Aguilar, 1997), 
in several works focused on different geological features 
of the north-central sector of the SMS (Fries 1960, 1966), 
and in more recent works that include radiometric age data 
(Meza-Figueroa et al., 2003; Moran-Zenteno et al., 2004). 
The present work is part of an ongoing study to examine 
the intrusive exposures, and represents the first study of the 
field, petrographic and geochemical characteristics of the 
plutons, with focus on the largest of the intrusive bodies, 
the Coxcatlan pluton. 

The Coxcatlan pluton, emplaced in Cretaceous lime- 
stone, is exposed over an area of about 18 km? and comprises 
at least two apparently coeval facies. Thin sections of 41 
samples from the Coxcatlan pluton were studied in a pe- 
trographic microscope in order to characterize mineralogi- 
cal and textural features. Modal analyses were carried out 
mainly to determine the amounts of quartz, K-feldspar and 
plagioclase present. Field observations indicate that most of 
the exposed pluton is a light gray colored, medium-grained 
rock speckled with variable amounts of dark minerals. In 
some parts of the pluton, crystals of plagioclase and alkali 
feldspar noticeably larger than the surrounding medium- 
grained minerals are common. These megacrysts are mostly 
<1 cm, but can reach up to 2 cm or more. Modal mineral 
analyses indicate that the Coxcatlan pluton plots mainly as 
monzogranite-granodiorite on a Streckeisen (1976) QAP 
diagram, with some samples extending into the quartz 
monzonite field (Figure 5); for simplicity, in this work 
the dominant facies of the Coxcatlan pluton is generally 
referred to below as a granodiorite. In this facies, biotite is 
ubiquitous (up to ~10%), minor hornblende (<2%) is often 
but not always present, and scarce pyroxene (up to 1%) and 
sometimes two pyroxenes were identified in a small number 
of samples. Titanite has not been observed. The granodiorite 
displays textural, mineralogical and slight geochemical 
variations. A finer-grained variety, exposed in the interior 
part of the pluton, is somewhat lighter in color compared 
to the dominant facies due to the smaller content of fer- 
romagnesian minerals (<5% biotite). Light colored felsic 
magmatic enclaves (see below) enclosed in the granodiorite 


could testify to an additional subsequent pulse of silicic 
magma. Xenoliths, rarely observed in the pluton, measure 
from 9 to 15 cm in length and have lithologies similar to 
those of the country rock. 

The second facies is represented by pyroxene-bearing 
magmatic enclaves, which are an important feature of the 
granodiorite. Mafic magmatic enclaves (MME) are ubiq- 
uitous and are finer-grained and darker than the host rock 
(Figure 4b, 4d, 4e, 4f, 4h); scarce felsic enclaves have been 
observed in the southern part of the pluton (Figure 4c). Most 
enclaves have ellipsoidal or rounded shapes, although some 
have irregular or lobate forms (Figure 4d); these forms sug- 
gest plastic deformation within a still fluid host. Contacts 
between the magmatic enclaves and the host rocks can be 
sharp with chilled margins or diffuse. Magmatic enclaves 
are usually <10 cm in length, but can reach up to 20 cm or 
more. Aligned belts of MME suggest dismembered dikes 
emplaced in a still fluid host rock (Figure 4e). Mafic mag- 
matic enclaves contain up to 10% quartz but typically less, 
plagioclase > potassic feldspar, and usually plot as quartz 
monzodiorites and monzodiorites in the QAP Streckeisen 
diagram (yellow circles in Figure 5). MME have a much 
higher content of ferromagnesian minerals with respect to 
the granodiorite. They contain abundant biotite (usually 
10-20%, reaching up to 25%) and pyroxene (usually <10%, 
but up to 15%); minor hornblende (generally <3%) can be 
present, but is observed as a reaction product of pyroxene 
in the transitional contact zone with the host granitoid. In a 
few enclaves elongate biotite grains show moderate to strong 
alignment, reflecting magmatic flow. Large plagioclase and 
alkali feldspar megacrysts are often observed in the MME 
and could have been transferred from the host granitoid 
when both magmas were still partially molten (Vernon, 
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Figure 5. QAP diagram for the plutonic rocks (Streckeisen, 1976). Red 
circles: Coxcatlan pluton; yellow circles: mafic magmatic enclaves hosted 
in the Coxcatlan pluton. 
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1986; Barbarin, 1990). Double magmatic enclaves (enclave 
within enclave) are frequently observed (Figure 4f). The 
presence of rounded to elongate enclaves, lobate contacts, 
double magmatic enclaves, and mafic and felsic enclaves 
side by side (Figure 4g) are features that suggest mingling 
(Vernon et al., 1988) in which pulses of magma of different 
compositions were injected into the host granodiorite. Mafic 
magmatic swarms are observed in the southern sector of the 
Coxcatlan pluton (Figure 4h). 

We have not observed any composite dikes in the 
area. Occasional thin aplite veins seven to eight cm wide 
are observed different parts of the Coxcatlan pluton, and 
a small aplite mass of approximately 40 cm is exposed in 
the northern sector of the pluton, along the highway to the 
northwest of the town of Coxcatlan. The aplitic rock is fine- 
grained and is constituted mainly by quartz and feldspar 
intergrowth in a graphic texture, with small amounts (~1%) 
of biotite and oxides. 

Recent widening of the highway at the northern 
edge of the Coxcatlan pluton has exposed small pegmatite 
lenses and a contact aureole between the granodiorite and 
Cretaceous limestone; yellow-brown garnet crystals in this 
zone measure 1—3 mm, but can reach 1.5 cm in diameter. 
X-ray diffraction in whole rock and clay-sized (<2 tm) 
fractions identified calcite, wollastonite, spessartite, vesu- 
vianite, and smectite in the contact aureole. XRD analyses 
were preformed with a Shimadzu XRD-6000, equipped with 
a monochromator, using Cu radiation. 

The Buenavista intrusion is a porphyritic subvol- 
canic body exposed in an area of approximately 2 x0.5 
km at the southern limits of the town of Buenavista de 
Cuéllar, and is limited to the east by ignimbrites and to the 
west by Mesozoic limestone, where it is associated with 
a Fe skarn deposit. It consists of phenocrysts of altered 
plagioclase crystals with disequilibrium features and em- 
bayed anhedral quartz grains, with scarce biotite in a fine 
granular silicified groundmass. Meza-Figueroa et al. (2003) 
report “°Ar/*°Ar plateau ages of 35.5 + 0.6 Ma and 34.7 
+ 0.6 Ma for two different experiments on biotite of the 
Buenavista intrusion. 

The Chautle pluton is a small intrusive body emplaced 
in Cretaceous limestone beds at a distance of ~3 km east of 
the caldera ring. It is a dark gray, fine- to medium-grained 
rock that contains plagioclase of variable size that includes 
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fine- to medium-grained crystals as well as large crystals 
showing disequilibrium features such as sieve texture, os- 
cillatory zoning, and resorbed margins. The Chautle pluton 
contains two pyroxenes (15—20%) with minor amounts 
of biotite (<5%), and has a dioritic composition, which is 
significantly more mafic than the dominant facies of the 
Coxcatlan pluton. 


GEOCHEMISTRY 


Major elements were determined by fused disk X-ray 
fluorescence spectrometry at the Laboratorio Universitario 
de Geoquimica Isotépica (LUGIS) of the Universidad 
Nacional Autonoma de México (UNAM), using a Siemens 
SRS-3000 spectrometer, following procedures outlined by 
Lozano and Bernal (2005). 

Whole rock major element analysis of volcanic units 
of the Tilzapotla caldera show a dacitic composition for 
the collapse-related Tilzapotla ignimbrite (~67-69 wt% 
SiO); post-collapse/pre-resurgence products have dacitic to 
rhyolitic compositions (~67—74 wt% SiO.), and lava flows 
related to resurgence (34.4 to 32.8 Ma) have andesitic to 
dacitic compositions (57-65 wt% SiO2) (Moran-Zenteno 
et al., 2004). Late stage volcanic rocks exposed between 
Buenavista and Coxcatlan were analyzed in this study; a 
hypabyssal unit (sample Cox15b) displays an andesitic 
composition with SiO, of ~57 wt% and a dacitic vitrophyre 
(sample Cox54) contains ~65% SiO, (Table 1). 

Geochemical analyses carried out on a group of 
intrusive rocks within the Coxcatlan-Tilzapotla elliptical 
structure display significant compositional variations in the 
late Eocene-early Oligocene plutonic record (Table 1). West 
of the caldera ring, the dominant facies of the Coxcatlan 
pluton ranges in composition from granite to granodiorite 
(~71 to 69 wt% SiO, (Cox26 and Cox7, respectively), but 
magmatic enclaves in the granitoids are generally much 
more mafic, and present a dioritic composition (SiO, = ~53 
wt%, Cox32). East of the caldera ring, the Chautle pluton 
(Chau1) is also more mafic than the dominant facies of the 
Coxcatlan pluton, with SiO, of ~55 wt%. 

In a geochemical study of the late Eocene-early 
Oligocene volcanic centers of the north-central SMS, 
Mori et al. (2012) suggest that the evolution of andesitic to 


Table 1. Whole rock major element analyses of magmatic rocks in the Coxcatlan-Tilzapotla area. 


Sample Rock Location SiO, TiO, 
Cox26 Coxcatlan pluton 99°27.98' 18°29.28' 70.97 0.32 
Cox7 —_Coxcatlan pluton 99°27.61' 18°28.30' 68.95 0.52 
Cox32. MME in Coxcatlan pluton 99°27.82' 18°29.78' 53.20 1.11 
Chaul  diorite 99°09.11' 18°25.95' 55.63 1.02 
Cox54__ vitrophyre 99°26.39' 18°28.99' 65.49 0.48 
Cox15b andesite 99°26.19' 18°28.63' 56.50 1.11 


ALO, Fe,.0; MnO MgO CaO Na,O K,O P,0O; LOI Total 
14.78 2.44 0.03 0.78 1.78 3.36 438 0.14 0.76 99.75 
15.15 4.15 0.07 1.27 2.87 3.34 3.52 0.14 0.35 100.35 
17.15 932 0.15 4.36 7.62 3.26 252 0.20 0.84 99.73 
17.31 823 0.13 5.13 742 2.88 155 0.19 0.43 99.91 
15.58 3.26 0.07 1.34 3.76 3.03 3.33 0.19 3.63 100.16 
18.47 7.28 O11 3.68 6.74 3.35 142 0.34 0.39 99.38 
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rhyolitic rocks associated with the Tilzapotla caldera was 
dominated by low-pressure fractional crystallization of 
pyroxene and plagioclase assemblages, with assimilation 
of Rb-enriched middle-upper crustal rocks to produce the 
dacitic and rhyolitic compositions. This interpretation was 
based on a coherent decrease of MgO, CaO, Al,O; and Sr 
concentrations with increasing silica contents, a negative 
correlation between Eu/Eu* and the differentiation index, 
and similar La/Yb and Sm/Yb ratios throughout the suite, 
which precludes the crystallization of garnet, a high pres- 
sure mineral with high partitions coefficients for heavy rare 
earth elements such as Sm and Yb. This is confirmed by 
high '*Nd/'“Nd ratios, low *’Sr/**Sr ratios, and a negative 
correlation between '*7Nd/“Nd and SiO . 


COMPARATIVE ANALYSIS OF THE 
GEOCHRONOLOGICAL DATA FROM THE 
PLUTONIC AND VOLCANIC ROCKS IN THE 
COXCATLAN-TILZAPOTLA REGION 


Geochronological dating was performed on distinct 
magmatic units not only to attempt to better constrain the 
evolution of magmas within the elliptical dome structure in 
the Coxcatlan-Tilzapotla caldera area, but also to explore 
pluton growth and the duration of the magmatic system prior 
to eruption. We present *’Ar/*°Ar mineral dating and U-Pb 


LA-ICPMS dating of zircon, in a number of instances from 
the same sample, of several intrusive bodies and extrusive 
units (Table 2). 


“Ar/*Ar geochronology 


Analytical methods 

Seven mineral separates were prepared from six 
whole rock samples collected in the Coxcatlan-Tilzapotla 
area. Four of the mineral separates were obtained from 
three samples collected in different parts of the Coxcatlan 
pluton: biotite-bearing samples Cox3 (biotite) and Cox10 
(biotite and alkali feldspar), and biotite-poor sample Cox26 
(alkali feldspar). Other samples dated are the Chautle diorite 
(plagioclase, sample Chau1), an intermediate volcanic rock 
(plagioclase, sample Cox15b), and the Atopula ignimbrite 
(sanidine, sample Cox20). 

The samples were prepared at the Mineral Separation 
Laboratory of the Departamento de Geologia at the Centro 
de Investigacion Cientifica y de Educacion Superior de 
Ensenada (CICESE), Ensenada, Baja California, Mexico. 
The basic sample preparation consisted of crushing and 
sieving, then rinsing with distilled water followed by 
98% acetone. The rock fragments were dried overnight at 
~60 °C. The minerals were then concentrated by Frantz 
magnetic separator to insure >99% purity. Samples were 


Table 2. Summary of *°Ar-*Ar and U-Pb geochronological data. 


Sample Location Mineral Unit Preferred Ar-Ar ages Apparent ?Pb/4U 
and K-Ar ages (Ma) zircon age (Ma) 

Cox3 99°28.03' biotite Coxcatlan pluton 36.28 + 0.23 ty 
18°30.12' 

Cox10 99°27.68' alkali feldspar Coxcatlan pluton 35.84 + 0.37 t, 35.75 + 0.28 
18°28.48' 

Cox10 99°27.68' biotite Coxcatlan pluton 35.76 + 0.13 t, 
18°28.48' 

Cox26 99°27.98' alkali feldspar Coxcatlan pluton 36.97 + 0.12 t, 36.01 + 0.36 
18°29.28' 

Chaul 99°09.11' —_ plagioclase Chautle diorite 34.54 + 0.55 ty 34.37 + 0.26 
18°25.95' 

Cox20 99°16.87' — sanidine Atopula ignimbrite 33.74 + 0.29 t, 33.18 + 0.28 
18°18.16' 

Cox15b 99°26.19' —_ plagioclase Hypabyssal andesite 30.77 + 0.90 t, 
18°28.63' 

Cox32 99°27.82' MME in Coxcatlan pluton 36.75 + 0.26 
18°29.78' 

Cox18 99°25.10' _ biotite* Buenavista intrusion I“ exp 35.5 + 0.6" t,* 35.17 + 0.18 
18°26.59' 2™ exp 34.7 + 0.6 

Cox19 99°10.64' —_ sanidine® Tilzapotla ignimbrite 34.26 + 0.09° one-step? 34.32 + 0.26 
18°22.62' 

Tz737 99°23.27' biotite‘ Tilzapotla ignimbrite 34.3 + 1.5° K-Ar® 34.23 + 0.20 
18°34.13' 


All age determinations from present work unless otherwise indicated. t,= plateau age; t,= isochron age; one-step = one-step 
fusion experiment on single crystals, a: Sample ABCO7 (Meza-Figueroa et al., 2003) collected from a location near Cox18; 
b: Sample Sol 2 (Moran-Zenteno et al., 2004) from same location as Cox19; c: Sample Tz145-01 (Moran-Zenteno et al., 2004) 


collected from a location near Tz737. 
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selected by handpicking under a binocular microscope 
from fractions ranging in size <710 to >350 um for samples 
Cox10 and Cox 20. Samples Cox3, Cox15b, Cox26 and 
Chau range in size from <350 to >250 um. As irradiation 
monitor, aliquots of sanidine FCT 2 (28.201 + 0.046 Ma; 
Kuiper et a/., 2008) were irradiated alongside the samples. 
The samples and monitor were irradiated with Cd-liner 
at the University of McMaster (Hamilton, Ontario) in the 
U-enriched research reactor receiving a dose of 30 MW. The 
4Ar-*°Ar experiments were conducted at the Geochronology 
Laboratory at CICESE using a Coherent Ar-ion Innova 
370 laser extraction system on-line with a VG5400 mass 
spectrometer. The irradiation monitors were fused in one 
step. All the samples were step-heated and additionally 
for sanidine Cox20 four one-step fusion experiments were 
performed. All the argon experiments were preceded by a 
blank measurement, with the five argon isotopes measured. 
Upon blank subtraction, the argon isotopes were corrected 
for radioactive decay, mass discrimination, calcium, potas- 
sium and chlorine neutron induced interference reactions. 
In processing the data, the decay constants and isotopic 
composition recommended by Steiger and Jager (1977) 
were used. Equations given by York et al. (2004) were used 
in all straight-line calculations. 


Results 

“Ar-*Ar results are summarized in Table 2 with full 
“Ar-*Ar data of all the experiments in the electronic ap- 
pendix Al, and details of each sample analyzed presented 
in appendix A2. Figures with age spectra and the **Ar/*Ar 
versus *Ar/“Ar correlation diagrams for each sample are 
presented in Figure 6. Note that all errors are reported at 
lo level. The errors in the integrated, plateau and isochron 
age include the uncertainty in the irradiation parameter J. 
Additionally, for the plateau and isochron ages, the good- 
ness of fit was included in the age uncertainty whenever the 
MSWD exceeded |. The integrated ages were calculated 
adding the fractions of the step-heating experiments. Plateau 
ages were calculated with the weighted mean of three or 
more consecutive fractions, which were in agreement within 
lo errors excluding the uncertainty in J. 

In general, the ““Ar-*°Ar and K-Ar age data are con- 
sistent with geologic relations and magmatic stratigraphy. 
Where available, biotite and alkali feldspar ages are concor- 
dant, suggesting that inherited-Ar is not a problem at this 
age range (e.g., Hora et al., 2010). The Coxcatlan pluton 
yields the oldest ages, ranging from 36.97 + 0.12 to 35.76 + 
0.13 Ma. At the 20 level these ages are still discordant sug- 
gesting a protracted cooling history in a body of significant 
volume. The Tilzapotla ignimbrite yields generally younger 
K-Ar and ““Ar-**Ar ages of 34.3 + 1.5 and 34.26 + 0.09 Ma, 
respectively, although the K-Ar age is relatively imprecise 
and is concordant within error with the younger Coxcatlan 
age. The ages of the subordinate units of the Chautle diorite, 
Buenavista intrusion, and Atopula ignimbrite are consistent 
with them representing post-caldera activity (Moran- 


Zenteno et al., 2004) and a general eastward younging of 
the magmatic system. 


U-Pb LA-ICPMS zircon geochronology 


Analytical methods 

The eight samples selected for zircon geochronol- 
ogy were crushed and standard techniques such as Wilfley 
separation and heavy liquids were used for mineral separa- 
tion. Zircon concentrates were observed under a binocular 
microscope and approximately 100 grains were selected 
from each sample, taking care to include all the morphologi- 
cal, shape and color variations. The selected grains were 
mounted in epoxy resin, ground by hand to expose half of 
the crystals, and imaged by cathodoluminescence using an 
ELM 3R luminoscope (Marshall, 1988). 

LA-ICP-MS U-Pb analyses were performed at the 
Laboratorio de Estudios Isotépicos (LEI), Centro de 
Geociencias, UNAM, using a Resolution M050 excimer 
laser workstation coupled with a Thermo Xseries quadru- 
pole ICPMS. The PleSovice standard zircon (ca. 337 Ma; 
Slama et al., 2008) was used as bracketing standard. All 
the analyses were performed employing a spot of 32 um in 
diameter and a fluence of ca. 6 J/cm’. A 1-sigma external 
reproducibility of 0.6% (7Pb/*8U ratio), 0.95% ?°’Pb/?°Pb 
ratio) and 1.1% (?°8Pb/?”Th ratio) was measured on such 
standard zircon, and quadratically propagated to the un- 
known measured ratios, to take into account the natural 
standard inhomogeneity. Further details of the analytical 
procedures are given in Solari et al. (2010). The obtained 
raw data were reduced, and errors propagated, using inter- 
nally developed software (Solari and Tanner, 2011). When 
needed, the common Pb contribution was corrected using the 
algebraic method of Andersen (2002), provided that the 7“Pb 
signal is generally swamped by the °“Hg interference and 
cannot be correctly evaluated by the current instrumentation. 
The obtained data were screened and filtered, following a 
discordance criteria (all the analyses which yielded more 
than 30% and less than -5% discordance were considered 
unreliable and thus discarded), and then an error criteria on 
isotopic ratios (all the analyses which corrected isotopic 
ratios yielded more than 5% | sigma error were discarded). 
Finally, further screening was performed observing rare 
earth element (REE) patterns, as well as some trace-element 
contents, such as P and La. The presence of analyses with 
high values (e.g., >1500 ppm P and >30 ppm La) strongly 
suggests that an apatite inclusion was hit during analysis, 
and were thus discarded. U and Th content were screened 
cautiously, because the higher their content, the greater the 
probability of crystal damage from radioactive decay, and 
thus of possible Pb loss. In due case, those analyses were 
also discarded. 

The concordia plots, deconvolution calculations, and 
age-error calculations were completed using Isoplot v. 3.7 
(Ludwig, 2008). Especially useful for the Paleogene ages 
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presented below are the mean *°Pb/**8U ages, given that 
the measurement of *°’Pb is problematic in those young 
zircons, and the consequent uncertainty in the °°’Pb/?°°Pb 
ratio is not a good indicator of geologically meaningful 
discordances. It is also common to observe scattering in the 
mean *°Pb/*38U ages, which yield MSWD (mean square of 
weighted deviates) values that are largely greater than one, 
an indication that the possibility of mixed age-populations 
exists. To try and recognize the different age components 
in those samples that showed an initial MSWD generally 
above 3, we used the deconvolution method implemented in 
Isoplot and based on the Sambridge and Compston (1994) 
mixture modeling method. Once the mixture components 
were recognized, their mean *°°Pb/*8U ages were plotted 
together with errors and recalculated MSWD in Figure 7. 
This figure includes U-Pb concordia plots (first column); 
deconvolution diagrams with frequency curves (red curves) 
or population density plots (PDP) and, for comparative 
purposes, *°Ar/*’Ar ages where available (second column); 
and mean age diagrams (on the right side of Figure 7). 


Results 

Details of analyzed zircons from the samples are given 
in appendices A3 and A4, and a summary of the apparent 
ages is given in Table 2. Details of the age data for each 
sample are in narrative form in appendix A2. Descriptions 
of most samples dated by this method are given above in 
the sections on “Geological setting and stratigraphy of the 
Coxcatlan-Tilzapotla area” and Appendix A2. 

Complex populations are revealed in population 
density plots or PDPs (Figure 7). Individual zircon ages 
in deconvolution plots (Figure 7, second column) reveal 
almost continuous distributions, but show slight breaks that 
suggest different pulses of zircon crystallization. This was 
confirmed by the application of the unmixing algorithm 
(Sambridge and Compston, 1994), by which two or three 
different populations with distinct weighted mean ages were 
obtained for each sample. 


Inherited ages and core-rim age variations in single 
zircons 

It is now well established that volcanic rocks often 
record a significant pre-magmatic history of inherited or 
xenocrystic “basement” zircons. Additionally, “antecrysts”’, 
zircons from an earlier confocal magmatic episode, are also 
common. The dominant population however, should be 


“autocrysts”, zircons that grew in the erupted magma (e.g., 
Miller et al., 2007, Folkes et al., 2011). The youngest zircon 
populations seen in the PDPs (Figure 7) correspond to these 
“autocrysts”, whereas the older populations correspond to 
“antecrysts”. 

Cathodoluminescence imaging revealed cores in 
some of the zircon grains, and oscillatory growth zon- 
ing was frequently observed. Cox 10 and 26 and Chaul 
have a significant number of grains ranging in age from 
Precambrian to middle Eocene that are consistent with the 
age of basement units in the region (e.g., Elias-Herrera et 
al., 2000; Talavera-Mendoza et al., 2005, 2007), or that can 
be explained by assimilation of sedimentary or metasedi- 
mentary crustal units. 

It is interesting to note that inherited zircon grains in 
the Coxcatlan-Tilzapotla samples are often of a different 
size and shape when compared to the main zircon popula- 
tion. The characteristics of the main zircon populations 
in the eight samples studied display certain variations 
(Table 3). In general, the late Eocene zircon autocrysts 
studied are colorless to pale yellow, have elongated pris- 
matic shapes, measure from ~130 to 340 um in length, 
and have aspect ratios of 2 to 4, with some exceptions; for 
example, many of the zircon grains are noticeably larger in 
Tilzapotla ignimbrite sample Tz737 (Table 3). In contrast, 
in Coxcatlan pluton sample Cox10 two small (~100 um), 
dark stubby inherited zircons yielded a slightly discordant 
age of 526 Ma and a concordant age of 65 Ma, and an 
unusually large, dark zircon has a core age of 167 Ma. 

The larger size of the cores and growth zones in some 
of the zircon grains allowed exploring age variations be- 
tween cores and rims. Analyses in two zircon grains of the 
Coxcatlan pluton yielded core-rim age pairs of 167 and 164 
Ma, and 92 and 35.6 Ma, respectively (Table A3, sample 
Cox10). Additional core-rim analyses in the samples yielded 
late Eocene age pairs that are interpreted as different pulses 
of the magma system that were injected into the magma 
chamber during different stages of development. Seven 
core-rim age pairs are indicated in the deconvolution plots 
of Figure 7, in which each pair has the same color error bar. 
Core-rim age variations in the seven pairs range from 1.6 
to 3.4 m.y. Core ages correspond to one of the earlier mean 
age groups obtained and are interpreted as antecrysts. Rim 
ages correspond to the younger mean age group defined 
by autocrysts. 

We turn our attention now to the autocrysts. These 


Figure 7 (following two pages). U-Pb zircon ages of samples dated by LA-ICPMS. In the first column on the left are U-Pb Tera-Wasserberg concordia 
diagrams, in the second column are deconvolution diagrams, and on the right are mean age diagrams. In the deconvolution diagrams, U-Pb zircon ages 
between 32 and 40 Ma are illustrated by diamonds with associated error bars; numbers on the vertical axis refer to the number of zircon analyses; red lines 
are frequency curves or population density plots (PDP), vertical black lines indicate the mean 7°Pb/**U ages with error (gray shaded areas), and the green 
dashed lines indicate *°Ar/*°Ar ages. Concordia diagrams, frequency curves, and errors were calculated using Isoplot 3.7 (Ludwig, 2008). Deconvolution 
diagrams were calculated using Isoplot, based on the unmixing model method developed by Sambridge and Compston (1994). Data-point error ellipses 
on concordia diagrams, and age error bars for mean ages on the deconvolution diagrams and mean age figures are 20 from the *°°Pb/*U weighted mean 
age. a) Coxcatlan pluton, Cox10; b) Coxcatlan pluton, Cox26; c) Mafic magmatic enclave, Cox32, embedded in the Coxcatlan pluton; d) Buenavista 
intrusion, Cox18; e) Chautle diorite, Chau1; f) Tilzapotla ignimbrite — Acapulco highway, Cox19; g) Tilzapotla ignimbrite — Cerro El Jumilar, Tz737; h) 


Atopula ignimbrite, Cox20. 
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record the most recent episode of magmatism and hold the 
keys to understanding the development of the magmatic 
system. 


DISCUSSION 
The record of zircon autocrysts 


Miller et al. (2007) regard the lower age trend of zircon 
arrays as autocrysts representing the best approximation of 
the crystallization age, as the last zircon growth or complete 
solidification of the pluton, assuming that Pb loss can be 
discounted. For the geochronological data obtained in the 
Coxcatlan-Tilzapotla area, the youngest mean ages in the 
deconvolution plots of Figure 7 are thus interpreted as the 
mean apparent age of the pluton or magma batch. The inter- 
pretation of the these ages as mean apparent magmatic ages 
is confirmed by the “°Ar/*°Ar mineral ages obtained, which 
are indistinguishable within uncertainty from the apparent 
age of the sample (Figure 7, Table 2). 

Weighted mean autocryst ages for the two samples 
of the Coxcatlan pluton are significantly older than those 
from the two samples of Tilzapotla ignimbrite. However, in 
detail the significance of this statistical discordance is offset 
by some interesting overlap in the individual. For instance, 
individual zircons from the Tilzapotla ignimbrite define a 
range in age from 36.6 + 0.4 to 33.3 + 0.5 Ma (lo), while 
Coxcatlan pluton zircons range in age from 39.5 + 0.6 to 
35.2 + 0.2 Ma. The autocryst populations of the two units 
show significant overlap in the PDPs, rim ages of Coxcatlan 
pluton autocrysts are younger or overlap with core ages of 
Tilzapotla zircons, and some individual Coxcatlan zircons 
that have ages concordant with Tilzapotla ignimbrite zir- 
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cons. We cannot reasonably attribute any of the overlap to 
Pb loss in the Coxcatlan pluton based on coherence of the 
data on the Tera- Wasserberg concordia plots. Thus, despite 
the statistical discordance of the weighted mean ages, 
we interpret these records to indicate significant overlap 
of zircon crystallization in the Coxcatlan pluton and the 
Tilzapotla ignimbrite. 

The Buenavista intrusion, on the west caldera rim 
and east of the Coxcatlan pluton, displays a wide range of 
zircon ages from 37.7 to 33.3 Ma. A major peak at ~35 Ma 
is defined in the PDP (Figure 7d), and yields a weighted 
mean age within error of the *’Ar/*’Ar biotite plateau ages of 
35.5 and 34.7 Ma obtained by Meza-Figueroa et al. (2003). 
Compared to the Coxcatlan pluton, zircon crystallization 
in the Buenavista intrusion began and ended later, but took 
place over a similar time period (~4 m.y.). 

In two samples (MME sample Cox32 and Buenavista 
intrusion sample Cox18), the youngest mean ages are 
represented by a small group of data with error bars that 
display little overlap with the main zircon population and 
are somewhat younger than the *°Ar/*’Ar ages (Figure 7). 
Again we cannot reject these as due to Pb loss as these 
appear coherent on the Tera-Wasserberg plots. They must 
represent real zircon crystallization ages and attest to dis- 
cordance between U-Pb zircon ages and “°Ar-*’Ar biotite 
ages in Cox 18 (no Ar-Ar age is available for Cox 32). We 
note that the “°Ar-*°Ar biotite age came from a different but 
closely located sample to Cox 18, that biotite can be open 
to Ar and K exchange during alteration, and that a skarn 
deposit is in this area. 

Ages for Cox 32 for the mafic enclave in the Coxcatlan 
pluton show significant overlap with the pluton data. The 
Buenavista intrusion sample Cox 18 also shows significant 
overlap. However, a significantly younger autocryst popu- 


Table 3. Samples dated by U-Pb zircon analysis: multiple analyses in single zircons and characteristics of analyzed grains. 


Sample Unit Core and rim ages 
(Ma) in single zircons 
Cox10 Coxcatlan pluton 167 and 164; 92 and 
35.6; 39.9 and 35.6 
Cox26 Coxcatlan pluton 
Cox32 Magmatic enclave in 38.5 and 36.6 
Coxcatlan pluton 
Cox18 Buenavista intrusion 36.8 and 33.7; 
35.1 and 33.5 
Chaul Chautle diorite 
Cox19 Tilzapotla 36.6 and 35.0 
ignimbrite 36.4 and 34.8 
Tx737 Tilzapotla 36.3 and 34.6 
ignimbrite 
Cox20 Atopula ignimbrite 36.4 and 33.0 


Zircon color Zircon size* Aspect ratio” 


(um) 

colorless to pale 200-250, 2 to 3.5, 
yellow up to 340 up to 4.5 
colorless 125-270, 2-3, 

up to 340 up to 3.7 
colorless to pale 165-280 2.6-3.4 
yellow 
colorless to pale 160-260, 2.5-4 
yellow up to 600 
colorless to very 130-260, 1.3-2.3 
pale yellow up to 400 
colorless to pale 130-330 2.5-4 
yellow 
colorless to pale 200-535 3-5.7 
yellow 
colorless to pale 150-260, 2.5-3.4 
yellow up to 400 


a: Average and maximum zircon size; b: Average and maximum zircon aspect ratio (length to width ratio). 


Caldera formation and progressive batholith construction, Coxcatlan-Tilzapotla area 263 


lation appears to be recorded in the Chautle diorite (37.0 
to 33.2 Ma) and the Atopula ignimbrite (37.6 to 32.8 Ma). 


Age patterns in the Coxcatlan-Tilzapotla area 


Significant overlap in all the zircon populations 
suggests continuous zircon crystallization in the Coxcatlan- 
Tilzpotla system as distinct batches of magma evolved 
over an approximately 6 m.y. period to build a significant 
volcano-plutonic complex. During this time, some of the 
complexity in the zircon populations could have been 
produced by recharge events, as evidenced by the Chautle, 
MME, and the young andesite. Hot magma intruding into 
a cooler, incompletely crystallized magma body preserved 
as supra-solidus crystal mush could result in changing 
the conditions of zircon saturation sufficiently to cause 
significant complexity. We are unable to resolve such events 
categorically but note the possibility of this for further 
investigation. 

A spatial trend of younger zircon ages is recorded 
to the east, not only in the Tilzapotla ignimbrite, which is 
the most voluminous volcanic unit and is associated with 
caldera collapse in the eastern sector of the dome, but 
also in the Chautle diorite and Atopula ignimbrite (Figure 
7). Ar-Ar geochronological data for these units are also 
younger. We interpret this to suggest a generally eastward 
migrating magmatic focus. It is of interest to note the U-Pb 
zircon apparent ages obtained for the Tilzapotla ignimbrite 
and Chautle diorite are the same within error, indicating a 
temporal relationship between the two units and suggesting 
a shared magmatic history (Figures 7e, 7f, 7g), in which the 
emplacement of mafic batches similar to and coeval with 
the Chautle diorite into the magma chamber could have 
triggered ignimbrite eruption and caldera collapse. 

Compared to the Tilzapotla ignimbrite, post-collapse 
ignimbrites (i.e., Atopula unit along the southern caldera 
limit) display a nearly continuous range in the zircon age 
data of ~4 m.y., starting at ~37 Ma and extending to a 
younger array at ~33 Ma, indicating the arrival of new 
batches of silicic magma into the chamber. *°Ar-*’Ar ages 
of the dacitic to andesitic lava flow units associated with 
the resurgence of the central block of the caldera floor and 
infill range from 33 to 32 Ma (Moran-Zenteno et al., 2004). 

The southeastward younging trend of the zircon 
ages reveals that batholith construction underneath the 
Coxcatlan-Tilzapotla dome progressed from west to east 
as pulses arrived in series into the upper crust. Individual 
pulses may have been emplaced sufficiently closely in time 
and space to have shared zircon crystallization histories, at 
least in part, as evidenced by overlap of rim and core ages 
in Coxcatlan pluton and Tilzapotla ignimbrite, respectively. 
The younging trend in the magmatic activity is illustrated in 
Figure 8. In the first stage, beginning at ~39.5 Ma, magma 
batches were emplaced in the western part and crystallized 
over a period of ~4 m.y. probably sustained by mafic 


Younging age trend in magmatic activity 


Vereen ee eso ee SE 


Coxcatlan 
39.5 — 38.0 Ma 


Vere SE 


Tilzapotla 
Coxcatlan ignimbrite  Chautle 
38.0 — 35.0 Ma 36.6 —33.3 Ma __ diorite 


34.4 Ma 


V—.ereoe?: eee SE 


Hypabyssal Resurgent 
Vitrophyre andesite events 
30.8 Ma 33-—32Ma 


Legend 
[| Vitrophyre 


LCL] Post-collapse dacitic to andesitic lava flows and dikes 
[| Tilzapotla ignimbrite 

LC] Coxcatlan pluton 

Cretaceous (Morelos/Huitzuco formations) 

[ Metavolcanic Taxco Viejo unit 


|| Mafic magma 
ES Magma chamber 


Figure 8. Sequence of sketches of the structural dome in the Coxcatlan- 
Tilzapotla area, illustrating the main magmatic episodes and the southeast- 
ward-younging age trend in magmatic activity. See text for explanation. 
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recharge. This produced an initial doming deformation in 
this sector. As the magmatic focus moved eastward, the 
most significant thermal anomaly developed from 36.6 
to 33.3 Ma, leading to the construction of a significant 
eruptible magma volume that extended the doming effect 
to the east. This doming in response to development of a 
major silicic magmatic system is consistent with classic 
geological observations of pre-caldera tumescence (e.g., 
Smith and Bailey, 1968) and recent thermomechanical 
models of caldera systems (Gregg et al., 2012). Continued 
intrusion of mafic magmas is evidenced by the Chautle 
diorite approximately concordant with eruption of the 
Tilzapotla ignimbrite and caldera collapse. Finally (~33-31 
Ma) aseries of silicic and intermediate injections of magma 
produced caldera resurgence and minor volcanic episodes 
outside the caldera. 


Batholith construction in the Coxcatlan-Tilzapotla 
area 


Recent advances of in situ U-Pb zircon dating of 
volcanic and plutonic rocks have shown that granitoid 
plutons (e.g., McNulty et al., 1996, 2000; Coleman et al. 
2004; Glazner et al., 2004; Matzel et al., 2006; Lipman, 
2007; Miller et al., 2007; Tappa et al., 2011) as well as the 
volcanic products of silicic magma systems (e.g., Schmitt 
et al., 2003; Reid, 2003; Bachmann et al., 2007a, 2007c; 
Folkes et al., 2011) record a multistage history of numerous 
pulses of magmas that accumulate over time and cannibal- 
ize their progenitors. We clearly see evidence for this in the 
Coxcatlan-Tilzapotla elliptical dome structure. We believe 
that the record here is one of incremental but continuous 
assembly of batholith from numerous pulses of magma. The 
nearly continuous age range observed for all samples sug- 
gest an approximately 6 m.y. period for the development of 
the system. This is consistent with, for example, Cretaceous 
batholiths in the western United States that have a duration 
of >5 m.y. for Mt. Stuart (Matzel et al., 2006) and >8 m.y. 
for Tuolumne (Coleman et al., 2004). 

The Coxcatlan pluton seems to represent an apophysis 
of a deep-seated plutonic complex, and displays textural 
and mineralogical diversity in addition to features that we 
interpret as key aspects in the understanding of the mag- 
matic system evolution. It is located in a segment of the 
structural dome where no significant volcanic conduits have 
been identified and where there is no caldera development. 
The overlapping age distributions (39.5 Ma to ~35 Ma) of 
zircons from different facies within the Coxacatlan pluton 
and observed mutual intrusion relationships documented 
between the mafic and grandiorite facies confirm that newer 
pulses were emplaced in an interval of ~4 m.y., without 
complete crystallization of previous pulses. 

The clear time-space relationship between the 
Coxcatlan pluton intrusion, the uplift in the Coxcatlan- 
Tilzapotla dome area, and the Tilzapotla caldera formation 


strongly suggest that these episodes are part of a continuous 
magmatic development driven and thermally sustained 
by mafic magma recharge. U-Pb zircon ages that span 
approximately 6 m.y., starting at ~39 Ma in the west part 
of the dome (Coxcatlan pluton) to ~33 Ma in the east 
(Tilzapotla and Atopula ignimbrites and Chautle pluton), 
as well as the compositional variations of pluton facies 
and volcanic units in the area can be taken as evidence 
that the dome uplift was the result of gradual construction 
of a batholith underneath due to the injection of different 
magma batches. 


What triggered the formation of the Tilzapotla 
caldera? 


The last decade has seen major developments in 
understanding of causes of intense episodes of explosive 
silicic volcanism, or ignimbrite flare-ups. Large storage 
regions of eruptible magma are characterized as crystal 
mushes that may contain over 50% of crystals and may 
themselves erupt or may be the source of residual liquids 
that are extracted and then erupted (Lindsay et al., 2001; 
Bachmann and Bergantz, 2008). According to numerical 
simulations of heat transfer in the upper crust, de Silva and 
Gosnold (2007) have shown that an elevated thermal flux is 
critical for the rheological conditions necessary to allow a 
large magma chamber to grow without erupting. These same 
conditions are necessary to maintain large melt-rich volumes 
of magma to feed caldera-forming eruptions (Annen, 2009). 
The recognition that such fluxes are most commonly found 
during magmatic flare-ups and are characteristic of batholith 
formation as well (e.g., Ducea, 2001) suggests that a view 
where caldera forming eruptions and pluton formation re- 
quire separate conditions is unwarranted. Large catastrophic 
caldera-forming eruptions that produce large ignimbrites are 
thus intimately associated with intrusive activity at depth 
(e.g., Smith, 1979; de Silva et al., 2006; Bachmann et al., 
2007b; Lipman, 2007). Extended plutonic activity would 
result in the development of crystal mushes that can lead to 
the accumulation of enough residual magma to produce a 
voluminous volcanic eruption. In such a scenario, the grow- 
ing “batholith” reaches a point where the “domed” roof is 
extended, faulted, and mechanically unstable. At a “tipping 
point”, the roof can be induced to collapse into the “cham- 
ber” producing a catastrophic caldera collapse (Gregg et al., 
2012) where most of the erupting crystal-mush accumulates 
as intra-caldera tuff (de Silva et al., 2006). The relative 
abundance of phenocrysts in juvenile pumice fragments 
as well as the occurrence of glomerocrysts with abundant 
plagioclase in the Tilzapotla ignimbrite suggest that the 
major volcanic eruptions related to caldera formation were 
associated with the disruption of crystal mushes. Compared 
to the plutonic units, the narrower age range of zircons in 
population density plots of the Tilzapotla ignimbrite (Figure 
7f, 7g) might be an indication of rapid development of the 
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eruptible magma in the central-eastern sector of the dome 
as the thermal conditions evolved (Figure 8). Overlapping 
zircon age ranges of the Tilzapotla ignimbrite and the small 
mafic Chautle pluton located near the caldera ring, as well 
as the Ar-Ar mineral age obtained for the Chautle pluton, 
which is within uncertainty of the zircon age, may suggest 
that injection of less evolved high-temperature magmas were 
instrumental in destabilizing crystal mushes (e.g., Bachmann 
et al., 2007b) and the roof and triggering caldera formation 
(Gregg et al., 2012). 


CONCLUSIONS 


Zircon U-Pb dating of late Eocene plutons and 
ignimbrites of the Tilzapotla caldera indicates that they 
represent the record a long-lived magmatic system 
genetically related to the construction of the large scale 
Coxcatlan-Tilzapotla volcano-plutonic complex and 
associated structural dome. As suggested by compositional 
variations and the multiple injections recognized in the 
Coxcatlan pluton, the dome is the superficial expression 
of a batholith that progressively grew from a series of 
injections of silicic magmas into the upper crust. “°Ar/??Ar 
mineral cooling ages obtained for the Coxcatlan pluton are 
generally concordant with the weighted mean *Pb/?*U 
zircon ages, suggesting shallow emplacement and relatively 
rapid crystallization. U-Pb age ranges, as well as mean 
ages of the Coxcatlan pluton with respect to ignimbrites 
of the Tilzapotla caldera and the Buenavista and Chautle 
plutons, indicate a southeastward migration of the thermal 
anomaly with the highest intensity being underneath where 
the caldera developed. In the southeastern part of the dome, 
catastrophic caldera collapse occurred due to the magmatic 
system reaching a mechanical threshold where the extended 
and domed roof above the growing batholith was unstable 
and eventually failed. Recharge with less differentiated 
high-temperature magma may have been the trigger for the 
eruption and caldera (roof) collapse as it disrupted crystal 
mushes and destabilized the magma chamber. 


APPENDIX A. SUPPLEMENTARY MATERIAL 


Appendices Al to A4 can be found at the journal 
web site <http://rmcg.unam.mx/>, in the table of contents 
of this issue. 
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